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Abstract. In this paper we study algorithms for tiling problems. We show that the 
conditions (Tl) and (T2) of Coven and Meyerowitz [5], conjectured to be necessary and 
sufficient for a finite set A to tile the integers, can be checked in time polynomial in 
diam(A). We also give heuristic algorithms to find all non-periodic tilings of a cyclic 
group Zjv- In particular we carry out a full classification of all non-periodic tilings of 
Z144. 
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1. Translational tiling 

Definition 1.1. We say that a subset A of an abelian group G (written additively) tiles 
G by translation when there exists a subset B C G such that every g £ G has a unique 
representation 

g = a + b, (a£A,beB). 

This situation is denoted by G = A(B B or G = B © A. The sets A and B will be called 
tiling complements of each other. 

Remark 1.1. In much of the bibliography, especially the earliest one and that with a more 
'— 1 algebraic focus, the tiling condition G = A® B is often called a factorization of the group 

In this paper we are primarily interested in the case when G is a finite group, typically 
ff~) a cyclic group Zjv = Z/(iVZ), or the group Z d , d > 1. 

m 

1.1. Periodicity. 

Definition 1.2. A set B C G is called periodic if there exists a nonzero g £ G such that 
B + g = B . Such a g is then called a period of B and clearly the set of all periods plus 
forms a subgroup of G. 

In case B C Z is a periodic set of integers, its least positive period is denoted by V(B). 

^ A periodic set B is clearly a union of congruence classes mod V(B) and we can write 

d B = B ®T{B)Z, 

where B C [P(B)Z] and we have used the notation 

[n] = {0,1,. ..,n-l}. 

It has long been known [20 j that if A is a finite set of integers and Z = A © B is a tiling 
of the integers then the set B C Z is necessarily periodic. A consequence of this is that 
the study of any such tiling with period, say, iV = V(B), reduces to the study of a tiling 
of the group Zjy, namely the tiling Z^r = A © B. 
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Non-periodic tilings of a group G = A © B are in some sense less structured and more 
interesting. Indeed, if h € G \ {0} is a period of B and H < G is the subgroup generated 
by h, then we can write B = B' © H for some B' C G/-ff (slight, harmless abuse of 
language here). The tiling condition G = A © B then becomes equivalent to the tiling 
condition G/H = A(B B'. This means that periodicity implies an immediate reduction to 
the complexity of the problem. Hajos [12] called a finite abelian group G a good group if 
in any tiling G = A © B one of the two sets A and B must be periodic. Sands |23| IM] 
completed the classification of all finite abelian groups into good and bad groups (for the 
entire list see also [TDJ Prop. 4.1]). The smallest cyclic group that is not good is Z72. 

1.2. Tiling in Fourier space. Cyclotomic polynomials. The condition G = A © B 
is clearly equivalent to 

(1) Yl Xa{9)xb{x - g) = 1, (xeG), 

g&G 

which we can rewrite as \A * XB = 1 ; with * denoting the convolution of two functions. 

Assume for simplicity that G is finite from now on. For a function / : G — > C its Fourier 
Transform [22] is defined as the function / : T — > C given by 

/(7) = E/(fb(5), 

T being the dual group of G. It consists of all characters 7 which are the group homomor- 
phisms G — > C. It is easy to see that / * g = f ■ g. Using this our tiling condition ([TJ can 
be rewritten as 

XA-XB = T = |G|x{ }- 
This is in turn equivalent to the obvious condition \A\ ■ \B\ = \G\ as well as 

(2) Z(£4)UZ(S) = G\{0}, 

where Z(f) = {x : f(x) = 0} is the zero set of /. So, to check that two sets A, B C G are 
tiling complements we have to first check the obvious condition \A\ ■ \B\ = \G\ and verify 
that the zero sets cover everything but 0. 

In the particular case where G = Zjv is a cyclic group, we have that T is isomorphic to 
G and the FT of / : G — > C now becomes the function on Z^r 

N-l 

f(k) = £ f(J)e-^' N , (k = 0, 1, . . . , N - 1). 

j=0 

Introducing the polynomial P(X) = ^2^=^ f(j)X^ it is often convenient to view the FT 
/ as the evaluation of P(X) at the N-th roots of unity 

-L,C ,0 ,...,0 

It is the roots of this polynomial P(X) that matter when one checks condition And in 
the case where / = \A is the indicator function of a set A C Zjv the polynomial P{X) has 
integer coefficients. The irreducible (over Z) factors of P(X) which are responsible for the 
zeros of P(X) at N-th. roots of unity are cyclotomic polynomials, which are defined as the 
minimal polynomials of roots of unity (see [5] for a brief introduction to the cyclotomic 
polynomials). 

A quick way to define the cyclotomic polynomials <F„(X) is by the decomposition of 
X N — 1 into irreducible factors 

X N -l = l[* d (X), 

d\N 

where the product extends over all divisors of N, 1 and N included. 
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The irreducible monic polynomial $ n (X) has as roots precisely the primitive n-th roots 
of unity, which are therefore algebraic conjugates of each other 

$ n (X) = Y\ {X-e 2mk/n ). 

0<k<n, (k,n)=l 

As a consequence, whenever an integer polynomial vanishes on a primitive n-th root of 
unity it vanishes on all of them. When vanishing of such a polynomial is the issue, 
therefore, the set of iV-th roots of unity is split into d(N) (= number of divisors of N) 
blocks, and the zero set of such polynomial P(X) among the roots of unity is a union of 
such cyclotomic blocks. 

Remark 1.2. Let us also observe that periodicity of a set A C Z^ can easily be detected 
on the Fourier side. Indeed, it is easy to see that, if N = ab, the set A C Z^y is a-periodic 
if and only if its Fourier transform xX{k) vanishes on all k £ {0, 1, . . . , N — 1} which are 
not multiples of b. 

1.3. The Coven- Meyerowitz conditions. Suppose A is a finite set of nonnegative in- 
tegers (we assume £ A) and write, as is customary, 

A{X) = Y,X a . 

Let Sa be the set of prime powers p a such that § p a(X) \ A{X). In [5 j Coven and 
Meyerowitz wrote down the following two conditions on a such a polynomial A{X). 

(Ti) A(l) = U 3e s A M^ 

(T2) If si, . . . , s m G Sa are powers of distinct primes then $ sl ... Sm (X) | A(X). 
They proved in [5 j that if (Ti) and (T2) hold for a set A then A tiles the integers by 
translation. Another way to say this is that there is a number TV" such that A tiles Ztv- In 
the converse direction they proved that (Ti) necessarily holds for any tile A. Regarding 
the necessity of (T2) for tiling it was proved in j5] that (T2) is also necessary for tiling 
when \A\ has at most two different prime factors. It was conjectured by Konyagin and 
Laba [H] that A is a tile of the integers if and only if both (T\) and (T2) hold. 

1.4. The computational status of tiling. The most basic computational problem of 
tiling is to decide whether a given finite set A in an abelian group G tiles the group. If 
one ignores, as a first approach, questions of complexity and restricts oneself to questions 
of decidability, one must assume G to be infinite (and, of course, discrete) for the problem 
to be meaningful. 

In a more general form of the problem, that of asking whether a given set of tiles can 
be moved around (by a group of motions) to tile M. d , tiling has long been shown to be 
undecidable. Berger |3] first showed this (it is undecidable to determine if a given finite 
set of polygons can tile M 2 using rigid motions). Many other models of tiling have been 
shown to undecidable. 

It is the case of tiling by translations and by a single tile that interests us here. 

When the group G is the group of integers Z the problem is decidable. This nontrivial 
fact follows from Newman's result [20] that every translational tiling of Z by a finite set is 
periodic. Although a bound for the period (2 diamj4 ) is given by Newman's theorem, this 
is not strictly necessary to deduce the mere existence of an algorithm to decide tiling by 
A (if, that is, we do not care about the running time of the algorithm). Indeed, given the 
finite A C Z, all we have to do is to start examining, by exhaustive search, whether A can 
tile the set {— N, . . . , TV} for larger and larger N. Tiling here means that we want to find 
a collection of non-overlapping translates of A which will cover {—N, ... , iV}. If the set A 
cannot tile Z then there is a finite N for which we will not be able to tile {— N, . . . , N}. 
This is a simple compactness (or diagonalization) argument. On the other hand, if A does 
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tile Z then, for large enough N, we will observe a tiling of {— N, . . . , N} which is "periodic" 
and which can, therefore, be extended indefinitely to the left and right to tile Z. 

The argument using periodicity given above is very general (see the introduction of (21] , 
where it is also proved that if a set admits a tiling of the plane with one period then it 
also admits a fully periodic tiling), works in all dimensions, and it is enough that periodic 
tilings exist. It is not necessary for all tilings to be periodic (which fails to be the case 
even in dimension 2 for as simple a tile as a rectangle). Thus, the so-called periodic tiling 
conjecture [18] (everything that tiles by translation can also tile periodically) implies 
decidability of translational tiling. This conjecture is still open in all dimensions d > 2. 

Already when the group is Z 2 the question of deciding if a given finite set i C Z 2 
can tile by translation is wide open, apart from the result of Szegedy [26] who gave an 
algorithm for the special cases of \A\ being a prime or 4. There are also algorithms for 
other special cases but these all have topological conditions [23 [TO] on the tile (e.g. to be 
simply connected). 

Let us now restrict ourselves to the decision problem of deciding whether a given A C "L N 
can tile Zjv by translation. We are interested to study the computational complexity of 
this problem regarding N as the parameter. In particular we'd like to have an algorithm 
which runs in time 0(N C ) for some fixed c > 0. Such an algorithm is still lacking though 
except when some arithmetic conditions on \A\ are imposed. In this paper we prove that 
if \A\ has at most two prime factors then we can decide if A is a tile of Z in polynomial 
time. This is so because we can decide the Coven-Meyerowitz conditions (T\) and (T2) 



(see £|1.3[) in polynomial time (Theorem 2.1 below). We also introduce a local version of 



the Coven-Meyerowitz conditions in ^3] which will allow us to give a similar algorithm, 
under the same arithmetic conditions, for when A tiles Zjv- 

On the other hand, some very similar problems to the decision of tiling (problem DIFF 
in [13]) have been shown to be NP-complete. This would suggest a lower bound in the 
computational complexity of the tiling decision problem. Were this decision problem 
to prove to be NP complete this would refute the equivalence of (T\ & T2) to tiling 
(conjectured in |14j). assuming of course P^NP. 

In the last section of the paper we turn to the problem of finding all non-periodic tilings 
of a cyclic group Z^r . As explained in subsection |1.1[ periodic tilings are less interesting 
because they can be considered as tilings of factor groups of Z^r. Finding many (or, indeed 
all) non-periodic tilings of Z^r could be a way of testing the Coven-Meyerowitz conditions 
and possibly producing counterexamples. Besides being mathematically interesting on its 
own right, this problem has a particular motivation in certain modern music compositions. 
The interaction of mathematical theory and musical background has been extensively 
studied in recent years [U 12) El [27] . 

H. Fripertinger listed out all non-periodic tilings of Z72 and Zios with the help of a 
computer search [HJ. In this paper we restrict our attention to N = 144, but we believe 
that the methods described here can be used to classify all non-periodic tilings for all 
N < 200. We also see from Section [4] that the number of non-periodic tilings grows at 
least exponentially with N, so that for large and highly-composite values of N this task 
gets hopeless. Also, the motivation for musical compositions does not extend beyond 
N = 200 — 300 due to obvious perceptional limitations. Let us also mention that in order 
to test the Coven-Meyerowitz conditions one will need to list non-periodic tilings of Z^r 
for values of N containing at least 3 different prime factors, such as N = 120, 180, 200, or 
even much higher. 

2. Deciding the Coven-Meyerowitz conditions 
In this section we describe in detail an algorithm which, given a set of integers 

AC {0,. ..,£>} 
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decides if A satisfies the (Ii) and (I 2 ) conditions of Coven and Meyerowitz [5]. 



Theorem 2.1. There is an algorithm to decide whether conditions (T\) and (T2) (see { 1.3) 
hold for a given A<^7L which runs in time polynomial in D = diam (A) . This algorithm 
therefore decides if A tiles Z when \ A\ contains at most two distinct prime factors. 

Proof. The algorithms consists of the steps given below. 



1. Compute all cyclotomic polynomials of degree up to D. 

This step is very easy to carry out in time polynomial in D using, for instance, the 
formula 

^n(x) = H(l-X n / d )^ d \ 
d\n 

where fj,(d) is the Mobius function. 

2. Determine the cyclotomic divisors of the polynomial A(x). 

Again, this is easily doable in time polynomial in D. Let Sa be the set of indices n = p^' , 
which are prime powers, such that & n (x)\A(x). Let pi,p 2 , ■ ■ -Pk be the different primes 
whose powers appear here, and let N{ > 1, i = 1, . . . , k, be the number of relevant powers 
of the prime pi . 

3. Test if condition (Ti) holds. 

Having computed the set Sa in the previous step, this amounts to checking the definition 
of condition (Ti). 

4. If (Ni + 1)(N 2 + 1) • • • (iVfc + 1) - 1 > D then answer that (I 2 ) fails. End. 

The explanation of this step is as follows. If (T2) is to hold then for any choice of prime 
powers (with respect to different primes) from the set Sa there is a different divisor of 
A(x), namely the cyclotomic polynomial & n (x) where n is the product of the chosen prime 
powers. Comparing degrees we obtain the inequality 

(JV 1 + l)(JV 2 + l)---(JV fc + l)-l< D. 



5. Check exhaustively if (T 2 ) holds and reply accordingly. 

This exhaustive search has to check (iVi + 1)(A^2 + 1) • • • (N/- + 1) — 1 possibilities for 
(T 2 ) to fail (select a power of each involved prime pi, with exponent from to Ni). This 
number is at most D (by the previous step) and therefore the total cost of this step of the 
algorithm is polynomial in D. 

The fact that this algorithm decides if A tiles Z when \A\ contains at most two distinct 
prime factors follows from Theorem B2 in [5 j which claims that, in that case, the conjuction 
Ti and T 2 is equivalent to tiling Z. □ 

3. The local Coven-Meyerowitz conditions 

In this section we give a local version of the Coven-Meyerowitz conditions, relevant to 
tiling a given cyclic group Z^v and not Z, and derive the corresponding properties they 
enjoy. 

Let iV be a positive integer. Define Sj[ be the set of prime powers p a \N such that 
Q p a(X) I A{X). Similarly define the conditions 

(If): A(l) = n^*.(l). 

(If): If G Sf are powers of distinct primes then $ sl ... Sm (X) | A(X). 
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In Theorem 3.1 below statements ([!]) and ^ correspond to Theorems A and Bl of [5]. 



Statement (|3j) means that the Coven-Meyerowitz conjecture is stronger than the corre- 
sponding local conjecture. Statement Q, corresponding to Theorem B2 in [5], means 
that we know the local condition is also true when at most two primes are involved 
in \A\. Finally, statement (|5j means that we have a polynomial time (in N) algorithm to 
decide the local Coven-Meyerowitz conditions and T% • 



Theorem 3.1. For any finite set ACZ and any positive integer N : 



(1) Tf and Tf A tiles Z 7 

(2) A tiles Z N ==> Tf 

(3) // (A tiles Z =>- Ti and T 2 ) then (A tiles Zjy =>• and T$). 

(4) If there are at most two distinct primes in \A\ and A tiles Z^r then and 
hold. 

(5) There is an algorithm which decides if and T 2 N hold in time polynomial in 
N. This algorithm therefore decides if A tiles Z^v when \A\ contains at most two 
distinct prime factors. 

Proof. The proof of ([TJ is essentially the same as that of Theorem A in [5] . When reading 
that proof it might help to observe that A(1)B(1) = lcm (Sa), or, in the local version, 
A(l)5(l) =lcm(Sf). 

For the proof of ^ we use Lemma 2.1 in (5j. According to that, if A tiles Z^r it follows 
that Sa contains only divisors of N. But A tiles Z^r implies that A tiles Z, hence, by 
Theorem Bl of 0, T-y is valid. But S A = S% in this case hence T^ is valid too. 

To prove (|3]), that the Coven-Meyerowitz conjecture implies the corresponding local con- 
jecture one first observes that if A tiles Z*n then it also tiles Z. By the Coven-Meyerowitz 
conjecture follows the validity of T\ and T 2 and, finally, observe that T 2 implies T 2 for all 
N. 

To show Q we have, if A tiles Z^r, that A tiles Z, hence T 2 holds (by Theorem B2 in 
[5]), which implies that TJ? holds. 

To prove ^ we merely repeat the algorithm of ^2] but replacing Sa by S^, etc. Because 
of Q this algorithm also decides if A tiles Z^v if \A\ contains at most two distinct prime 
factors. □ 



4. The number of non-periodic tiling complements 

One might expect that the number of non-periodic tilings of Z^v is small. However this 
is not the case. 

Theorem 4.1. There are arbitrarily large N and non-periodic tilings Z^r = A® B, such 
that there are additional distinct non-periodic tiling complements B\ , . . . , of A, with 
k > e c ^ , C a constant. 

Proof. The following sketch of a proof relies on a construction given in [T7] where it was 
used to prove that there are tilings of the integers with a tile contained in {0, . . . , D} but 
of period at least CD 2 . 

We take iV = 2 • 3 • 5 • p ■ q where p and q are two different large primes, roughly of the 
same size ~ \/N. We view the group Zjy as 

Z N = Z 3p x Z 5q x Z 2 

as shown in Fig. [T] We take the set A to be the "3 x 5 rectangle at the origin" . The 
set A can then tile Z^r in a non-periodic way by first tiling the lower and the upper layer 
in the ordinary way and then perturbing a row of the lower layer and a column of the 
upper layer, as shown in Fig. [TJ Each such perturbation manages to destroy one of the 
two periods that existed in the ordinary tiling (for a detailed proof see [T7] ) . 
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Figure 1 . The group Z 3p x Z 5q x Z 2 



If one wants to create many different non-periodic tilings of Zjv with A one does as 
follows. 

(1) Keep the upper layer the same (with one perturbed row only). 

(2) Keep half the lower layer the same: rows indexed from 5q/2 to 5q are left unper- 
turbed (refer to Fig. [T]) . 

(3) Perturb the remaining rows of the lower layer arbitrarily, subject only to the re- 
striction that at least one of them is perturbed by a non-zero amount. 

It is clear that the number of such complements is at least 

Each such tiling is non-periodic: suppose g = (a, b, c) S "L^ p x x 7L<i is a period. The 
two layers cannot be swapped by a translation by g, as one has perturbed rows but the 
other one does not, so c = 0. The upper layer is mapped to itself by a ^-translation. This 
can only happen if b = as the perturbed column is unique. Finally a = since the lower 
layer cannot move in the direction of 7L^ q and be mapped into itself: the perturbed row 
closest (in the "direction of motion") to the unperturbed block of ~ 5q/2 rows cannot 
move to a perturbed row. 

Even if one considers two tilings which differ by a translation identical the number would 
still remain exponential in N. □ 



5. Finding all non-periodic tilings in a cyclic group 

As explained in Section [T] there is a particular motivation coming from music compo- 
sitions to list out all non-periodic tilings of a cyclic group Ztv, for relatively small values 
of N. H. Fripertinger [8] achieved this task for N = 72,108. In this section we settle 
iV = 144, and we believe that the algorithms described here are likely to work for other 
values like iV = 120, 180,200,216. On the other hand we see from Theorem 4.1 that this 
task becomes hopeless for large and "highly composite" N. 
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We cannot offer any rigorous mathematical statements here with regard to the com- 
plexity of the algorithms used. We provide heuristics, observe that methods seem to work 
remarkably well and we are able to achieve a full classification for N = 144. 

Remark 5.1. Throughout this section we rely on results of [5j. Our algorithm works 
only when tiling implies the TJ? condition (see ^jj. In particular it works when either N 
has at most two distinct prime factors, such as N = 144 = 2 4 • 3 2 , or has three distinct 
prime factors but one of them appears to the first power, such as N = 120 = 2 3 • 3 • 5 or 
N = 180 = 2 2 • 3 2 • 5. The reason for this is that in any tiling Z^r = A © B in such a case 
either \A\ or \B\ has at most two distinct primes factors. 

A basic ingredient of our algorithm is the "fill-out procedure" performing the following 
task: given a set A in a group G, list out all tiling complements B of A. We cannot give 
any precise statement on how efficiently this task can be achieved except to note that it 
performs really well in practice. However, let us mention here that the method described 
below was already successfully used in (6] in refuting a conjecture of Lagarias and Szabo on 
a necessary condition for the existence of universal spectrum, and it could well be useful 
in tackling other question related to tiling. The credit for this heuristic algorithm goes to 
P. Mora in [B] . 

The Fill-out Procedure. 

Assume A C Z^r is given (in fact, the procedure works for any finite abelian group G, 
but we restrict our attention to cyclic groups). The task is to find all tiling complements 
B of A. As usual, we make the normalization assumption that £ A, B. 

We will build up, by adding elements one by one, all tiling complements B of A. For 
this, we move in the space of packing complements P of A, i.e. sets P = {p\, . . . p m } 
such that all translated copies A + pj are disjoint from each other. We are trying to grow 
the sets P so that they become tiling complements. Our exploration starts with the set 
P = {0}. At any given step we will attempt to extend P by one new element in a particular 
manner (which is described in detail in the next paragraph). If no such element exists we 
backtrack. We describe our algorithm "explore" below, and note that explore(P) finds 
all tiling complements B containing P. We invoke it as explore({0}) with the list of tiling 
complements of A being initially empty. The recursive procedure explore is described in 
Fig.[2} 

The efficiency of this procedure is due to the heuristic of expanding a given packing 
complement P of A by adding to it a copy of A which covers an element x that is the 
most restrictive in the sense that there are few possible ways to cover it. This is achieved 
by ranking the yet uncovered elements x in increasing order of the function r{x) = r(P, x) 
that counts in how many ways an element x can be covered by adding one more copy of 
A to the given packing P. Clearly, if r{x) = for some x then x cannot be covered in any 
admissible way, so that P can definitely not be extended to a tiling complement (it is still 
possible that P could be extended by several further elements to form a larger packing 
complement P' , but it can never grow to become a tiling complement due to x never being 
covered). It also often happens that r(x) = 1, i.e. we are forced to place the next copy of 
A at a certain position in order to cover x, and once this copy is placed, there is even less 
room in Z^r, and it is likely that there will be another x such that r(x) = 1, or at least 
small. 

With this fill-out procedure at hand we can now describe the algorithm to list all non- 
periodic tilings of Z144. We merely describe the algorithm here, as the full documentation 
(or listing out the numerous arising sets) would be far too lengthy for this paper. The 
results and documentation are available online at |16j . 
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explore(P) 

(1) If P has already been explored, return. 

(2) Mark P as explored. 

(3) If P is a tiling complement of A, add it to the list of tiling complements and 
return. 

(4) Compute the function r : TL^ \(A(B P) — > N, defined by r(x) = number of ways to 
add an extra copy of A to A © P so that x is now covered and the new copy of A 
does not intersect A © P. 

(5) Rank all x S \ {A © P) according to r(x): 

r(x\) < r(x 2 ) < < r(xk), 

in such a way that if r(xi) = r(xj+i) then Xi < Xi+i (for this purpose of ordering 
the elements of Z^v are thought of as {0, 1, ... N — 1}). Here k = \1*n \ {A © P)|. 

(6) If r{x\) = 0, return. 

(7) Consider x\. 

• Let A + J/1, . . . , A + y r t xi ) be the copies of A that can be added to A® P in a 
non-overlapping way and that contain x\. 

• explore (PU{yt}), . . ., explore (P U {y x ( ri )})- 

Figure 2. The explore recursive procedure 



Normalizing conditions. We always assume that € A, B. Also, it is clear that in 
terms of tiling questions a set A = {0, a%, ct2, . . . } is equivalent to its own translated copies 
A — d, A — a>2, ... all containing zero. It is customary to include only one representative 
from these equivalent copies, and we are going to follow this tradition (so that the numbers 
described in the last section correspond to this normalization). 

The Algorithm. 

Step 1. The prime powers dividing 144 are 2, 4, 8, 16, 3 and 9. In any tiling A®B = Z144 
the cyclotomic polynomials corresponding to these prime powers must divide exactly one 
of A(x) and B(x), according to condition T\ of [5 . Therefore we first make a list of all 
possible partitions {H,H C } of the elements {2,4,8,16,3,9}. There are 32 such partitions 
(note that {H,H C } and {H C ,H} are the same). Our task is to decide for each partition 
whether any non-periodic tilings correspond to it. 

Step 2. Certain partitions produce only periodic tilings due to condition T2 of [5]. 
Indeed, as an example, take the partition {{2,4}, {8, 16, 3, 9}}. Due to condition T2 of [5] 
the cyclotomic polynomials &2i(x), $72(2;), $4s(x), $144(2;) must also divide B(x), together 
with &s(x), §>ig(x), $3(2;), $9(2;). From this it follows that the support of xb is contained 
in the subgroup {0, 2,4,..., 142} and this automatically makes B periodic. Therefore, in 
this step, discard all partitions which imply automatic periodicity of either A or B by 
condition T2. 

The remaining partitions are: 
{{3,4,8}, {2,9,16}}, {{3,4,9}, {2,8,16}}, {{3,4,16}, {2,8,9}}, {{3,8,9}, {2,4,16}}, 
{{3,8,16}, {2,4,9}}, {{4,8,9}, {2,3,16}}, {{4,8,16}, {2,3,9}}, {{4,9,16}, {2,3,8}}, 
{{8,9,16}, {2,3,4}}, {{2,3,4,8}, {9,16}}, {{2,3,4,16}, {8,9}}, {{2,3,8,9}, {4,16}}, 
{{2,3,8,16}, {4,9}}, {{2,4,8,16}, {3,9}}, {{2,4,9,16}, {3,8}}, {{2,8,9,16}, {3,4}}, 
{{3,4,8,16}, {2,9}}, {{4,8,9,16}, {2,3}}, {{2,3,4,8,16}, {9}}, {{3,4,8,9}, {2,16}}, 
{{2,3,4,9}, {8,16}}, {{2,3,4,8,9}, {16}}. 
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We remark here that most of these partitions will not produce non-periodic tilings. 
If there were further theoretical considerations upon which certain partitions could be 
discarded, then it would probably be possible to apply our algorithm for higher values of 
N. For instance, it could be true that if either \A\ or \B\ is a prime power then the tiling 
must be periodic. This would allow us to discard several further partitions. 

Assume we are dealing with a fixed partition P = {H,H C }. The natural approach is 
the following: by the structure theory given in Section 4 of [5] one can list out all subsets 
A C Z144 such that A tiles Z144 and $h( x ) divides A(x) for all h G H. 

We remark here that besides using the structure theory of [5] there is an alternative 
way of listing out these sets A. Namely, take the least common multiple L of the prime 
powers appearing in H. By the Remark following the proof of Theorem A in [5] one can 
construct a universal tiling complement B in Zj, of any tile C in Z^ such that &h(%) 
divides C{x) for all h £ H. The set B is constructed by taking the product of certain 
cyclotomic polynomials as described in the proof of Theorem A of [5] . Next we apply our 
fill-out procedure to B in Z^, and hence list out all possible sets C in Z^. Finally, all the 
desired tiles A in Zjy must reduce to one of these sets C modulo L. Therefore, to list out 
the sets A we must shift the elements of sets C by multiples of L in all possible ways. 

Then these sets A need to be grouped into equivalence classes according to the zero set 
of the Fourier transform of xa- A similar procedure must be done for subsets B. Then 
one needs to select the equivalence classes which correspond to non-periodic sets A and 
B (this can be read off from the zero-sets of the Fourier transform, or directly from any 
particular representative of the class). Finally, if a class of non-periodic A's and a class 
of non-periodic i?'s are such that the union of the zero-sets of their Fourier transform 
contains Z144 \ {0}, then any representatives of these classes tile with each other in a 
non-periodic way. 

For some partitions this approach can indeed be carried out. However, the problem is 
that for certain partitions the number of tiles A and/or B provided by the structure theory 
of [5j can be overwhelmingly large. Take the example {{2, 3, 8, 16}, {4, 9}}. The structure 
theory of jS] provides two different tiles A±,A 2 of cardinality 24 modulo 1cm (2, 3, 8, 16) = 
48, such that A(x) is divisible by each of "£{2,3,8,16} ( x )- The tiles A modulo 144 are therefore 
obtained as sets that reduce to A\ or A2 modulo 48. There are 2 • 3 23 such sets. 

Step 3. To overcome this problem, we only list one type of the tiles (either ^4's or B's, 
whichever are fewer). In the above example we list the tiles of type B. The structure theory 
of |5 provides six different tiles B\,...Bq of cardinality 6 modulo 1cm (4, 9) = 36, such 
that B{x) is divisible by each of ^{4,9} (x). The tiles B modulo 144 are therefore obtained 
as sets that reduce to any of Bi,...B$ modulo 36. There are 6 • 4 5 = 6144 such sets. 
We can easily sort them into equivalence classes according to the zero-set of their Fourier 
transforms. There are five classes C%, . . . , C5 corresponding to non-periodic tiles B. They 
are characterized by the cyclotomic divisors $>h(x), with h being in C\ = {72,36, 18,9,4}, 
C 2 = {36,12,9,4}, C 3 = {36,18,9,4}, C 4 = {36,18,12,9,4}, C 5 = {72,36,18,12,9,4}. 

Step 4- Next we discard those equivalence classes Cj which automatically make any 
tiling complement A periodic. Consider, for example, C\ above so that our tile B has cy- 
clotomic divisors ^{72,36,18,9,4}- Then any tiling complement A of B must have cyclotomic 
divisors ^{2,3,6,8,12,16,24,48,144} (and possibly others), which automatically make A periodic. 

Step 5. For any remaining equivalence class Cj we take a representative Bj, and run 
our "fill-out procedure" to find all tiling complements A of Bj. The point is that Bj is 
non-periodic, so that heuristically we expect not too many tiling complements of Bj to 
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exist (although this heuristic breaks down if N is large and highly composite, as shown 
by Theorem 4.1). Once we have all tiling complements A of Bj we can select the non- 
periodic ones, if any. Let Nj denote the collection of all non-periodic tiling complements 
A of Bj. Then we conclude that all sets A in Nj tile with all representatives B of the class 
Cj. We repeat this procedure for all j and we arrive at a list of all non-periodic tilings 
corresponding to the fixed partition P. In the specific example above, Step 4 already 
discards all classes Cj so that Step 5 becomes unnecessary. 

We repeat Steps 3-5 for all partitions P listed in Step 2, and we arrive at a list of all 
non-periodic tilings of Z144. 

As a last normalization step, we must recall that a tile A = {0, 01, . . . } and its shifted 
copies A — ax, etc. are considered equivalent, and keep only one representative for each 
tile. And, naturally, the same applies to the tiles of type B. 



Step 6. Exceptional cases. We have been able to execute the steps above for all 
partitions listed in Step 2 except for {{2, 4, 8, 16}, {3, 9}}. In that case the number of tiles 
provided by the structure theory of |5] is too large for both A and B, and therefore the 
numerical search would take several days. We choose to deal with this case by invoking 
some results from the literature. We will show that either A or B must be periodic. 

Note that \A\ = 16 and \B\ = 9 by condition T\ of [Hj. Recall the following. 

Theorem 5.2. (Sands |25J If A © B = %n and N has at most two prime factors p,q 
then either A or B is contained in the subgroup p7L^ or qL^. 

Assume \ A\ is contained in a subgroup. Then the subgroup must be H3 = {0, 3,6,... 144} 
due to obvious cardinality reasons. Then B is a union of three parts Bq, B±, B2 according 
to residues modulo 3. All three parts must have cardinality 3, and A © Bq gives a tiling 
of the subgroup H3. Note that H3 is isomorphic to Z48. If A is not periodic then Bq 
must be periodic (all tilings of Z48 are periodic). Therefore we must have Bq = {0, 48, 96}. 
Similarly, B\ and B2 must be periodic by 48, and therefore so is B. 

Assume now that B is contained in a subgroup. Then the subgroup must be H2 = 
{0, 2,4, ... , 144}. Then A is a union of Aq and A±, the even and odd parts of A, and 
Aq © B = i?2- Note that H2 is isomorphic to Z72. If B is non-periodic then Aq must be 
periodic (the only non-periodic tilings of Z72 contain tiles of cardinality 6 and 12, while 
\B\ = 9 and \Aq\ = 8). In fact, Aq must be periodic by 72, because the periodicity 
subgroup contained in Aq must divide the cardinality \Aq\ = 8. By the same reasoning we 
obtain that A\ is also periodic by 72, and therefore so is A. 



Summary of the results for Z144. 

For the partition {{2, 8, 9, 16}, {3, 4}} there is a set T\ containing 36 non-periodic tiles of 
cardinality 24 (type A), and a set T2 containing 6 non-periodic tiles of cardinality 6 (type 
B). Each tile in T\ tiles with each in T2. The 36 tiles in T\ fall into two different equiv- 
alence classes corresponding to the following cyclotomic factors: ${144,72.24,18,16,9,8,2} > an d 
${144,72,18,16,9,8,2}- The 6 tiles in T2 fall into one equivalence class corresponding to the 
cyclotomic factors ${48,36,24,12,6,4,3}- We present here one representative for T\ and T2 
(one can then easily recover the full sets by applying the "fill-out procedure" to these 
representatives) : 

• {0, 17, 20, 23, 28, 29, 40, 48, 53, 59, 65, 68, 76, 88, 89, 95, 96, 101, 116, 124, 125, 
131, 136, 137} is in T\ and 

• {0, 32, 58, 90, 112, 122} is in T 2 . 
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For the partition {{4, 9, 16}, {2, 3, 8}} there is a set Ti containing 6 non-periodic tiles 
of cardinality 12 (type A), and a set T 2 containing 324 non-periodic tiles of cardinal- 
ity 12 (type B). Each tile in T\ tiles with each in T 2 . The 6 tiles in T\ fall into one 
equivalence class corresponding to the following cyclotomic factors: ${144,36,18,16,9,4}- The 
324 tiles in T 2 fall into two equivalence classes corresponding to the cyclotomic factors 
${72,48,24,18,12,8,6,3,2}, and ${72,48,24,12,8,6,3,2}- We present here one representative for T x 
and T 2 (one can then easily recover the full sets by applying the "fill-out procedure" to 
these representatives): 

• {0, 34, 40, 46, 48, 58, 88, 96, 106, 118, 130, 136} is in T x and 

• {0, 16, 29, 44, 57, 73, 80, 93, 108, 109, 124, 137} is in T 2 . 

For the partition {{2, 4, 9, 16}, {3, 8}} there is a set T\ containing 8640 non-periodic tiles 
of cardinality 24 (type A), and a set T2 containing 3 non-periodic tiles of cardinality 6 
(type B). Each tile in T\ tiles with each in T 2 . The 8640 tiles in T\ fall into three different 
equivalence classes corresponding to the following cyclotomic factors: ${144,36,18,16,12,9,4,2}, 
${144,36,18,16,9,6,4,2}, an d ${144,36,18,16,9,4,2}- The 3 tiles in T 2 fall into one equivalence class 
corresponding to the cyclotomic factors ${72,48,24,12,8,6,3} We present here one represen- 
tative for T\ and T 2 (one can then easily recover the full sets by applying the "fill-out 
procedure" to these representatives): 

• { 0, 9, 17, 26, 27, 34, 39, 40, 48, 51, 57, 65, 74, 82, 88, 96, 99, 105, 111, 113, 122, 
123, 130, 136} is in Ti and 

• {0, 36, 64, 80, 100, 116} is in T 2 . 

The partition {{3, 4, 8}, {2, 9, 16}} is the most interesting. There are sets T\ and S\ con- 
taining respectively 6 and 156 non-periodic tiles of cardinality 12 (type A), corresponding 
to the cyclotomic factors ${72,48,36,24,12,8,6,4,3} and ${72,48,24,12,8,6,4,3}, respectively Also, 
there are sets T 2 and S 2 containing respectively 12 and 48 non-periodic tiles of cardinality 
12 (type B), corresponding to the cyclotomic factors ${144,36,18,16,9,2} an d ${144,18,16,9,2}, 
respectively. Sets in T\ tile with all sets in both T 2 and 52- Sets in S\ only tile with sets 
in T 2 . We give a representative for each of T\, S\,T 2 , S 2 : 

• {0, 18, 28, 44, 54, 64, 80, 82, 98, 108, 118, 134} is in T lt 

• {0, 16, 30, 44, 58, 74, 80, 94, 108, 110, 124, 138} is in Si, 

• {0, 33, 40, 45, 48, 57, 88, 96, 105, 117, 129, 136} is in T 2 , and 

• {0, 27, 30, 35, 60, 72, 75, 83, 102, 123, 131, 132} is in S 2 . 

The other partitions listed in Step 2 do not produce non-periodic tilings. 
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